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We show that the stress-energy tensor has additional terms with respect to the ideal form in 
states of global thermodynamic equilibrinm in fiat spacetime with non-vanishing acceleration and 
vorticity. These corrections are of quantum origin and their leading terms are second order in 
the gradients of the thermodynamic fields. Their relevant coefficients can be expressed in terms 
of correlators of the stress-energy tensor operator and the generators of the Lorentz gronp. With 
respect to previous assessments, we find that there are more second order coefficients and that 
all thermodynamic functions including energy density receive acceleration and vorticity dependent 
corrections. Notably, also the relation between p and p, that is the eqnation of state, is affected 
by acceleration and vorticity. We have calculated the corrections for a free real scalar field - both 
massive and massless - and we have found that they increase, particularly for a massive held, at very 
high acceleration and vorticity and very low temperature. Finally, these non-ideal terms depend on 
the explicit form of the stress-energy operator, implying that different stress-energy tensor of the 
scalar held - canonical or improved - are thermodynamically inequivalent. 


I. INTRODUCTION 

It is common wisdom that the form of the relativistic stress-energy tensor in a thermodynamic equilibrium state 
has the ideal form: 


= {p + p)uf^u^ - pg^^'' 

where p and p are the energy density and pressure, thermodynamic functions of temperature T and chemical potential 
p, and u a constant four-velocity. In quantum statistical mechanics, the above expression corresponds to the renor¬ 
malized mean value □ of the quantum stress-energy tensor operator, built from local quantum fields, with the density 
operator: 

p=(l/Z)exp[-/3-P + CQ] (1) 

where /? = (1 /T)m is a constant inverse temperature four-vector (or, simply, four-temperature), with T = 
being the proper (or comoving) temperature and u the constant four-velocity, () = p/T is the ratio between proper 
chemical potential and proper temperature, P the four-momentum operator and Q an internal conserved charge: 

= tr(^'"'^(a:))ren = 7tr(f'"'^(a;) exp[-/3 • P-bCQ])ren = {pFp)u^^u'' - P9^'' (2) 

The form ([^ is dictated by the symmetries of the density operator Q which is traslationally invariant and isotropic 
in the rest frame where (3 = (1/P)(1, 0 ). 

However, the density operator is not the only form of global thermodynamic equilibrium, which is, in general, 
a state where the entropy S = —tr(plogp) is constant. For instance, it is well known [1] d] that in non-relativistic 
quantum mechanics the operator: 

p= (1/Z) exp[-P/ro + w7/To] (3) 

where Tq is a constant global temperature □ H the hamiltonian and Jz the angular momentum operator along some 
axis z, represents a globally equilibrated spinning fluid with angular velocity uj. Similarly (see sect. 0 the operator: 

p = (1/Z) exp[-P/To + aKz/To] (4) 

Kz being the generator of a Lorentz boost along the z axis, represents a relativistic fluid with constant comoving 
acceleration along the z direction and it is still an equilibrium distribution. These two cases belong to a more general 


^ For free quantum fields, by renormalization we mean the use of normal ordering in the stress-energey tensor operator 
^ The global temperature To is a temperature measured by a thermometer at rest with the external observer. In general, it differs from 
the proper temperature T measured by a comoving thermometer 
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class of thermodynamic equilibria which, in special relativity, are characterized by a four-temperature j3{x) field 
fulfilling the equation 


d^l3„ + = 0 


(5) 


which means that the four-temperature is a Killing vector field. We will show how for such thermodynamic equilibrium 
states, with the appropriate treatment in quantum relativistic statistical mechanics, the ideal form of the stress-energy 
tensor gets quantum corrections - vanishing in the /i —>■ 0 limit - whose leading terms are proportional to the squared 
gradients of /3, which in turn can be expressed in terms of the acceleration and the vorticity fields (see sect. IV 
for definitions): 

= {p + p)uf^u'' - ( 0 ( 0 ^) -t 0(uj^) + 0{auj)) 


As we will see, these corrections are normally tiny but they can become relevant under specific circumstances and, 
moreover, they are not microscopic in the sense of being relevant only at very small scales. 

The appearance of these terms is somehow in contrast to the widespread belief that deviations from the ideal form 
([^ can only arise in presence of dissipative processes. In fact, the existence of such terms has been pointed out by a 
classification of second order gradient corrections of the stress-energy tensor in conformal hydrodynamics Hi] also by 
means of kinetic theory and some coefficients, denoted as thermodynamic in view of their survival at equilibrium, 
have been calculated in ref. [6] for conformal field theories. 

In this paper, we show that the occurrence of non-dissipative corrections to the ideal form of the stress-energy 
tensor is a general fact which is related to the very notion of equilibrium in quantum relativistic statistical mechanics. 
Moreover, these corrections result from the expansion of the density operator and their form is not assumed a priori 
like in the Landau-frame based gradient expansion. At equilibrium, they have simple and suggestive expressions as 
correlators of the stress-energy tensor with the generators of the Lorentz group. The proper energy density expression 
is also modified, as well as the relation between energy density and pressure, that is the equation of state. It is an 
almost straigthforward consequence that these corrections will extend to a curved spacetime. 

The paper is organized as follows: in section|ll]we obtain the form of the density operator of general thermodynamic 
equilibrium in quantum statistical mechanics in flat spacetime. In section in we discuss the relation between local 
observables and the local value of the four-temperature field. In section [TV] we derive the form of the corrections to the 
ideal form of the stress-energy tensor as a perturbative expansion. In section|V]we calculate those quantum corrections 
in free scalar field theory. Finally, in sections |VI| and |VII| we discuss the most important physical consequences and 
draw the conclusions. 


Notation 

In this paper we use the natural units, with h = c = K = 1. 

The Minkowskian metric tensor is diag(I, —1, —I, —1); for the Levi-Civita symbol we use the convention 
We will use the relativistic notation with repeated indices assumed to be summed over, however contractions of indices 
will be sometimes denoted with dots, e.g. u-T ■ u = Operators in Hilbert space will be denoted by a large 

upper hat, e.g. T while unit vectors with a small upper hat, e.g. v. The stress-energy tensor is assumed to be 
symmetric with an associated vanishing spin tensor. 


II. EQUILIBRIUM IN RELATIVISTIC QUANTUM STATISTICAL MECHANICS 


A general covariant form of the density operator in relativistic quantum statistical mechanics extending the eq. Q 
was first proposed, to our knowledge, in refs. mil]: 


P = ^ exp 


- f dE^ 


( 6 ) 


where S is a spacelike 3D hypersurface. This form can be obtained maximizing the total entropy with the constraints of 
given energy-momentum and charge densities at some specific ’’time” of the hypersurface E, see the detailed discussions 
in ref. [8| and more recently in refs. [Ulini- The density operator (j^ is therefore especially suitable to describe local 
thermodynamic equilibrium — that is a situation where the thermodynamic parameters temperature, velocity field 
and chemical potential are a function of space and time — in a quantum relativistic framework. The operator (§ 
will not maintain its form under the unitary time evolution and cannot thus represent the actual quantum state 
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in the Heisenberg representation. However, it is time independent or, equivalently, independent of the integration 
hypersurface S if the divergence of the integrand vanishes and in this case the is the density operator of a 
thermodynamic equilibrium state. For conserved stress-energy tensor and current this condition leads to the request 
m that C is a constant and /3 a Killing vector field fulfilling eq. (§ (with partial derivatives replaced by covariant 
derivatives if necessary). 

The density operator ([^ is also well suited to describe thermodynamic equilibrium in a general curved spacetime 
possessing a timelike Killing vector field. In Minkowski spacetime, which we will be dealing with in this work, the 
general solution of the eq. ([^ is: 

/?- = 6" + (7) 

where & is a constant four-vector and zu a constant antisymmetric tensor, which, because of Q can be written as an 
exterior derivative of the (3 field: 




Hence, the general equilibrium form in flat spacetime of the density operator <© reads: 


P = ^ exp 


-b^P^ + + CQ 


( 8 ) 

(9) 


where the J’s are the generators of the Lorentz transformations: 


= 



Therefore, besides the chemical potentials, the most general equilibrium density operator in Minkowski spacetime can 
be written as a linear combinations of the 10 generators of the Poincare group with 10 constant coefficients. The most 
widely known case is the one with 13 = b and w = 0, that is eq. 0. what we define as homogeneous thermodynamic 
equilibrium. The rotating global equilibrium in eq. ([^ can be obtained as a special case of eq. 0 setting: 

bp = (1/To, 0, 0, 0) = {u:/To){gifj,g 2 v - giug2p) 

where w has the meaning of a costant angular velocity [1]. Similarly, the form Q can be obtained by setting: 

bp = (1/To, 0,0,0) ^p,^ = (a/To)(50^531/ - 53 m5oi/) 


In the latter case, the contravariant components of (3 read: 

( 3 ^ = ^ (1 + az,0,0,at) 


( 10 ) 


thus the unit vector j3 is the velocity field of a fluid with constant comoving acceleration along the field lines (for the 
field line going through z = 0, the comoving acceleration is a). 


III. MEAN VALUES OF LOCAL OPERATORS 


Suppose we want to calculate the mean value of a local operator 0{x) (in the Heisenberg picture) with the density 
operator 0: 


0{x) = {0{x)) = tiipO{x))ren = ( exp 


- f dS^ - C?) 


0{x) 


( 11 ) 


If is a general field, there is no compelling reason why, at a given point x, the mean value 0{x) should be simply 
equal to the same value at the homogeneous global thermodynamic equilibrium with an uniform four-temperature (3 
equal to its value in the point x, that is I3{x). For instance, the stress-energy tensor in the point x does not need to 
be of the ideal form 0 with u = /3(x) and p = p{l3'^,() p = p(/3^, O if 13 is not constant. In fact, its tensor structure 
in 0 is determined by the symmetries of the density operator Q, which is obtained from 0 provided that P is 
constant. 
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Nevertheless, one can imagine that if /3 and C are sufficiently slowly varying in space and time, 0{x) will be mostly 
determined by the values of the fields /3 and C around the point x [9]. More specifically, the distance over which the 
thermodynamic fields like /3 vary should be much larger than the typical thermal correlation length, which is governed 
by the microscopic parameters of the theory and the temperature itself. This can be shown by recasting the fields in 
the integrand of the eq. (111 as /3 = /3(x) + <5/3 and ( = ({x) + 6 C,, so as to obtain: 


0 {x) = ^tr 

1 


exp 


= —tr exp 


-P^ix) / dS^ T'"'" + C(a:) / 02.^ j 
Jn JT. 

-P.{x)P'' + ax)Q - dE^ - 5C?) 




0 {x) 


Hence, applying the linear response theory to the exponent above: 

0{x)^{d{x))p(,)- ["dz [ dE^(2/) ({d{x)f>^^{y + tzl3{x)))hi.)-{d{x))pi.){f>^''{y 


■iz(3{x)))i3(^^)) 6/3, 


-m(2/) { {Oix)j>^{y + izP{x))) 0 (^) - (O/x))■ 




( 12 ) 


+ dz dE,. 

Jo Jy. 

where x and y both lie on the hypersurface E. The symbol ( stands for the (renormalized) mean value calculated 
with the homogeneous equilibrium density operator in eq. (fll). Particularly, the ( stands for the mean value 
calculated with a fixed four-temperature (and C) equal to the value of the /3 (and C) fields in the point x. The 
formula (12) just expresses the aforementioned concept, namely that the local equilibrium value of the operator 0 {x) 


is determined by the local values of the thermodynamic fields with corrections depending on quantum-statistical 
correlations between operators in different points. These correlations - hence the integrand function in eq. ( |12[ ) - are 
significant over microscopic distances I dictated by the mass, temperature and coupling constants of the theory, which 
are supposedly much smaller than the macroscopic distance L over which 6/3 and 6 (/ appreciably vary. This condition 
is usually referred to as hydrodynamical regime and in this regime terms beyond the linear in (5/3 and 6C in eq. (12) 
contribute less and less as they are expected to be suppressed with higher powers of l/L. 

Under these circumstances, it is possible to expand the thermodynamic fields in the eq. 0 into a Taylor series 
about the point x. For a general /3 field, this method makes it possible to find an approximate expression of the local 
thermodynamic equilibrium operator (§ as a function of /3{x) and its derivatives [5]. For the special case of global 
equilibrium, with constant Q and /3 a Killing vector field 0 . one can recast the operator 0 so as to have in the 
exponent the value of the four-temperature in the point x: 


0 {x) = -tr 

1 


exp 


exp 


-h^P^ + + CQ 


0 {x) ) = -tr 


exp 




M)pt^ 


2^' 


T/ijy 


<Q 


0 {x) 


-/3^(x)p^ + Jr + CQ 


0(a 


(13) 


(14) 


where we have used the angular momentum operators around the point x: 

T(a:) = exp[ix • P] being the translation operator. 

The calculation of mean values (13) is the main purpose of this paper, and, specifically, when O = We will 
consider the term in w as small compared with the terms involving /3 and (/ and expand accordingly. Thus, the leading 
term in the above equation will be singly the homogeneous equilibrium one wit h fo ur-temperature equal to its value 

with u = /3{x). We will see in the sect. 


in the x point, that is the expression 


IV 


that the lowest order corrections 
to the ideal form are of the second order in vc and that they are of either quantum or quantum-relativistic nature as 
they vanish for /i —>■ 0 or /i/c —>■ 0. 

Note that /3{x) is required to be a future-oriented timelike vector in order to get a finite value for most observables 
at the lowest order of the /3 expansion. This condition cannot be fulfilled everywhere for the expression ^ ii zu ^ 0. 
For instance, for the rotating global equilibrium (|^, it is easy to check that: 

P = X x) 

which becomes spacelike when ||a;k x x|| > 1, that is when the velocity exceeds the speed of light. Similarly, for the 
operator (|^, the /3 field (10) is future-oriented timelike only in the Rindler wedge defined by the light cone of the 


point (0,0, 0, —a/T). Therefore, the validity of our calculations will be limited to the physical regions where the /3 
field is timelike and with positive time component even though the operator 0 written with the constants b and vj 
does not make this limitation apparent. 
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IV. PERTURBATIVE EXPANSION FOR THE STRESS-ENERGY TENSOR 


The goal of this section is to provide an expansion in w for the mean value of the stress-energy tensor in the general 
form of thermodynamic equilibrium: 


(15) 

/ ren 

Indeed, w is an adimensional tensor in natural units and it the has, in general, very small components. To understand 
its physical meaning, it is very useful to decompose it into two spacelike vector helds, each having three independent 
components, projecting along a timelike vector. A physically interesting choice is u = /3 = in the regions 

where /3 given by eq. Q is timelike. We can then decompose w as follows: 

+ a^u'^ - a''u^ (16) 


= -tr ( exp 


+ + CQ 


where, by definition: 


a^{x) = w^{x) = (17) 

Note that a and ic, unlike w, are not constant and they are both orthogonal to rt, hence spacelike. The physical 
meaning of a and w vectors can be shown starting from the eq. ([^. Because of Q and ([^ at equilibrium one has: 

fiiy — 


whence: 


We can now take the scalar product with and conclude that: 

u,d''y/^=D^ = 0 

which tells us that, as expected, at the thermodynamic equilibrium the comoving temperature along the flow lines 
does not change and 9^/?^ = V^/3^, where: 


^ fj, — 


Thereby, the a vector simply becomes: 


= - 


T 


(18) 


that is the acceleration field divided by the proper temperature. Note also that, being 8^(3^ = 0, one has: 


0 = u''{d^l3p + = ap+ pT 


2V^ 


T 


2^2 


(19) 


Likewise, it can be shown that w corresponds to an angular velocity divided by a temperature, for, by using 




= ^u,d,up = ,Up = 


T 


( 20 ) 


being uUp, as it is known in literature, the local vorticity vector. Restoring the physical constants, 

one then has the adimensional four-vectors: 


http 

cKT 


^p 

KT 


( 21 ) 


These numbers are, for the vast majority of physical systems, much less than 1 and a perturbative expansion in w 
of the eq. (15) is then feasible. They can give rise to relevant corrections if the implied additional terms to the ideal 
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stress-energy tensor are some sizeable fraction thereof or when these terms are comparable to the viscous tensor. 
According to the eq. (21), this happens at for very large accelerations or very low temperatures. Hence, let us define: 


TZ(w) = exp 


-/3^(x)p^ + -tn^,jr+cQ 


= exp 




exp[CQ] 


( 22 ) 


where, in the last equality, advantage has been taken of the supposed commutation of the charge operator Q with 
both the P’s and J’s. At the second order in vu one can write: 


7l{zu) = -I- + o(t77^) 

and, by using the Poincare group commutation relations, it can be shown that (see Appendix A): 
•^( 0 ) _ q-P-p+CQ 

= \{e-P-P+iQ ^ 

j^{2)fj,upa _ + ^q-P-P+CQ jjPjjPp'^p'^ _ ^q-P-P+CQ pfj-g’^pp'^ 

16 8 12 


(23) 


(24) 


where the curly bracket expression { , } stands for the anticommutator. 

By using the eqs. (23) and (24), the mean value (15) can be expressed as an expansion in zu with coefficients which 
are calculated at the homogeneous thermodynamic equilibrium: 

tT{K{w}} ^ 


pa '^\t 


iRe(Jr Jr;f^^(x))^(,) + \l3»{x)l3\x){P'^pP,fP''{x))p^,) - ^/3^(a:)ff"'^(P-f'^'^(x))^(,) 


--Re(Jr;T^"(x))^(,)(Jr)/ 3 (.) 


+ o(w^) 


(25) 


where we have used the relations for two hermitian operators A, B: 


Re{AB)=^-{{A,B}) 


tlm{AB) = -{[A,B]) 


and the notation has been introduced: 


{A-B) = {AB)-{A){B) 


for the correlator between A and B. The terms in eq. (25) containing P can be readily calculated taking the derivative 
of (T )/3 with respect to /3. Indeed: 


/3V" = -P‘^g^^{PP,fnp 


d(3a dPr 


({p^pp,fp'')p - {p^-,fnp{pip - {p^\fp'')p{p-)p) 


(26) 


Note that the two rightmost terms in the last equation vanish once multiplied by Wp^vJxr for, being (P)/? oc /3, they 
contain the symmetric combination /S'''/?’’ or jSPfi'^. 

All mean values in eq. (25) involving angular momentum operators can be rewritten in a form which makes it 
apparent that their dependence on x is only through the value of the four-temperature, by taking advantage of the 
translational invariance of the density operator. For instance: 


and similarly for the others, where eq. (14) has been used. Then, it is convenient to decompose the tensor J into two 
spacelike vector operators the same fashion as for w in eq. (16): 




( 27 ) 
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being 


= UpJP>^ 




The operators and are simply the generators of the rotation and boosts with respect to the reference frame 
with time direction u. Using the invariance by rotation (in the hyperplane orthogonal to u), parity and time reversal, 
which are assumed to hold for our hamiltonian, one readily obtains that (see Appendix B): 

Re(jP-;f/^-(0))^(^^ = 0 = 0 


Therefore, plugging the decomposition ((27| into the eq. (25), and using the relations 0 and ( [2^ and after the 
removal of the vanishing terms, the eq. ([25| can be written as: 


TP'^ix) = (T^"(x))^(,) + \apa„Re{KPk'^-fP''{0))p^,^ + ]^WpW,Re{JPJ'^ {0)) 


1 


92 


1 


d 


The derivative terms are easy to work out by using ([^ ; they will give rise to expressions involving the thermodynamic 
functions pressure, energy density and their derivatives, that is specific heats. On the other hand, the correlators 
in eq. (28) cannot be expressed in terms of known thermodynamic functions. In fact, they can be written as linear 


combinations of new thermodynamic coefficients which can be expressed in turn as correlators of specific components 
of the stress-energy tensor and angular momentum or boost operators J and K, that is 

-Re({.A^, k'^)] = -u^u'^AP'^ktiT, () + AP^AP’^kg^T, Q + {A^'^AP'' + A'^'^App - ‘^AP'^AP'^)k,{T, C) 

2 3 

iRe({>, J'^}; f'^‘^( 0 ))^(,,) = -uPu''AP'^jt{T, () + AP''AP’^jg{T, () + -k A'^'^A^^ - ^A^^A^'")j«(r, C) 

z o 


Re({if^ J"};T'^'^(0));3(,) = [upU^e^P'^'' + u''u,e^P'^P)l4T,0 
where 

AP'' = gP’' - uPu‘ 

is the projector onto the hyperplane orthogonal to u and 

3 

kt{T,C)=Re{k^k^;f°kO))T 


(29) 


jt(r,C)=Re(j3j3;f™(0))T 


HTX) = \Re^{k^k^;k\0))T 

1 ^ 

j,(T,C) = -Re^(j3j5;f"(0))r 


fc3(T,C)=Re(K'if2;fi2(0))r 

js{T,C)=Re{kj^;kkO))T 


UT,0=Re{{K\k}-T^kO))T 


(30) 


In the eq. ( [30| ) the notation ( )t has been introduced meaning that the expectation value is calculated in the rest 
frame where (3 = (1/T, 0 ). The derivation of eqs. (29) and (30) can be found in Appendix B. 

Finally, after having worked out the derivatives of the stress-energy tensor and using the eqs. (29) and § in the 
eq. (28), one obtains: 


TP''{x) = {p — a^Ua — vkUu,)uPu'' — {p — — vkD^)AP'' -|- AaPa'' + WwPw" -I- G{uPj^ + ^Pu") + o(ti72) (31) 


where p,p are the usual homogeneous thermodynamic equilibrium functions energy density and pressure, and the 
functions U, D, A, lU, G read: 


r. 1 - 1 ■ 

— 2^^ ~ 


Uw = ^jt Da = ^{p + p) + ^kg - ^k, 

^{p+p) + ks W =js 


( 32 ) 










The vector 7 ^^ in eq. (311 is defined as: 


= (a • = e^^'P^w^apU, 


(33) 


where the vj decomposition (16) has been used, 


As it can be seen from eq. (31), the stress-energy tensor has corrections to its ideal form which depend on quadratic 
combinations of the two vector fields, a and w arising from the decomposition of the exterior derivative of the 
temperature four-vector j3. At thermodynamic equilibrium, according to the previous discussion and the eqs. ( 211 , 


they are proportional to the acceleration and angular velocity (or vorticity) , so that the eq. (31) can be rewritten 
in the most suggestive fashion by restoring the natural constants as: 






(MYu 

\KT) “ 


P + 


/i|a| 

cKT 


Dn 






+.4 1411 


A'"" 


(34) 


where |a| = and |w| = y/—ujpU}>^ and a, uj are the corresponding unit vectors. In the expression (34) 

the adimensional scales Ha/cKT and huj/KT involving acceleration and vorticity have been separated from the 
thermodynamic functions U, A,W, D,G having the same dimension as p and/or p making it easier to appreciate the 
size of the correction to the ideal form. 


A. Relation with other second-order hydrodynamical coefficient calculations 


The appearance of extra terms in the stress-energy tensor at thermodynamic equilibrium with respect to its ideal 
form has, needless to say, several physical consequences. The presence of non-dissipative quadratic corrections in 
the vorticity and gradients of temperature (hence accelerations at equilibrium, according to eq. was pointed 

out in ref. [Sill] and, being non dissipative in nature, defined as thermodynamic in ref. |6|. The calculation of such 
coefficients has attracted much attention lately (see m for a recent review), especially in conformal held theories 
[IMS] with different techniques jHIITj (see also ref. [IH])- The coefficients that we have denoted by Da, D^, A and 
W are in the following relation with those known as ^ 3 , ^ 4 , A 3 , A 4 in literature: 


A 

rp2 

rp2 



W 

rp2 


— ^3 


i/a 

rp2 


(3A4 - 9^4) 


(35) 


Remarkably, the number of coefficients quoted in ( |34[ ) is larger than envisaged in ref. [3111 and the reason is that 
we did not assume, as it is usually done in the Landau frame, that the proper energy density p has the same 
functional dependence on the temperature as at homogeneous thermodynamic equilibrium. This assumption proves 
to be incorrect, and the extra coefficients cannot be reabsorbed by a redehnition of temperature, as it will be discussed 
and shown in sect. Ell 

Before tackling these issues, it is necessary to calculate the coefficients U, D, A, W, G in some instance and we will 
do it for the simplest case of a real scalar free field. As it will be clear from the calculations shown in the next section, 
they all have a classical expression in the massive case, and, as a consequence, all the corrections in eq. (34) to the 
ideal form turn out to be of quantum origin, as they vanish in the —>■ 0 limit. 


V. THE FREE SCALAR FIELD 


The goal of this section is to calculate the coefficients in eq. (34) for a free real scalar field. This implies C = 0 
in the density operator 0> yet it is quite easy to extend the obtained results to the charged case with 0 in the 
Boltzmann limit of distinguishable particles. 

The theory is described by the Lagrangian density: 
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By adding the super-potential 




a class of stress-energy tensors can be obtained as Noether currents associated to space-time translations. Although 
they are explicitely dependent on the parameter they differ from each other by a divergence: 




(36) 


thus they lead to the same generators of the Poincare group. For ^ = 0 the tensor is the so-called canonical stress- 
energy tensor, while for ^ = 1/6 the tensor is the so-called improved stress-energy tensor |19j . In the translationally 
invariant homogeneous equilibrium Q all mean values of local operators are independent of x, thus the divergence 
in the above expression vanishes, hence p and p do not depend on In fact, as we will show, this is not true in the 
case of generalized equilibrium and the correlators in eq. (30) are explicitely dependent on 


At the very beginning, it should be pointed out that in principle one should use normal ordering in the calculation 
of the mean values of T in a free field theory to subtract zero point infinity. However, this is not needed in the 
calculation of a correlator such as (T“^(0))r because :T:=T — (0|T|0) (T being a quadratic operator in the 
fields) so that the vacuum term cancels out in the subtraction (J'’‘^T“^( 0 ))t — (J^'^)t(F“^( 0 ))t. 

The basic tool we need in order to carry out the calculation is the free field n-points Wightman thermal function: 


W^\xi,X2,...,Xn) = {-ipixi)'ijjix2) ... '0(a:„))r 


which can be written in terms of 2-points thermal functions according to a version of the Wick theorem 
for thermal field theory. For an even n: 


suitable 




(n) 


n 




, X j _ 1, Xj _|_ X, . 


.) 


i =2 


while, if n is odd, W^\xi ,..., Xn) = 0. In the case of a free real scalar field, the 2-points Wightman thermal function 


reads: 


W^\x,y) = j d'^ke [ 6 »(fc°)-f nT(|A:°|)] - m^). 


where 


nrie) 


1 

qC/T _ ^ 


is the Bose-Einstein distribution. We then define 

y, z) = {f^^^x) fp^iy) f “'5 (z))t - (fP’^ix) fp%y))T {f^^{z))T , 

which can be calculated with the point-split procedure as 

q-P^\P'^\°‘l^[x,y,z) = 0 “^ yVx\xi,X 2 ,yi,y 2 ,zi,z 2 ) - ^(xi,2:2, j/i,2/2) ^(2^1,^2) 

where 

©r = {(1 - - 2idP^d:^ + \gP'' [(4C - l)d,p + m 2 ] . 

The general expression of the correlators is then: 

(>‘'>‘";f“'^(0))T = J d^xd^y[xPyPT^'''°'"'°'^(x,y,0) - x''yPTx^'°‘"''"^{x,y,0) 

-xPy'^T^''''°^''°‘'^{x, y, 0 ) -f {x, y, 0 )] 

with = ?/° = 0 because the J’s are time-independent. 


( 37 ) 
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Out of the 15 different diagrams stemming from the contractions of the 6 -point Wightman thermal function, in 
T some are cancelled by the subtraction term, leaving only the 12 diagrams in which T°‘^[z) is not a disconnected 


component. Since {J^^)t = 0 , in the eq. (37) the remaining 4 disconnected graphs in T do not contribute to the 
result. Therefore in (37) we can replace T with its connected subset of 8 diagrams and we get: 

1 


C. 


with 


/ii/|p(7|a/3 


(x,y,0) = 


(27r)9 


J d^k d4pd^ge-*(*^+P)^ e-i(g-k)y q) Sik"^ - rr?) 5 {p^ - m^) 5 {q^ - rr?) 

X [ 9 {k°) + nT{\k°\)] [ 0 (/)+nT(|/|)] [ 0 {q°)+nT{\q°\)] , 


V^''\p'^\°‘^ik,p, q)={-{l- 2£,){k^^p'' + p^k'') + 2^{k>^k'' + p^p'') - g^'' [(4^ -l)k-p + 2^{k^ + p^) - m^] } 

X {(1 - 2C){kPq'^ + qPk"^) + 2^(Ffc'" -b gV) - g"" [-(4C -l)k-q + 2^{k^ + q^) - m^] ) 

X {-(1 - 20{p"qP + gV) + 2e(pV + - 5“^ [(4^ - 1)P ■ <7 + " "i'] } • 


The thermodynamic correlators in eq. (30) can be found by selecting the suitable indices in eq. (37). For instance, 
for the kt correlator: 


kt{T) = / d^xd^y 40100100 


^O—yO—Q 


(38) 


Using 


and 


J d^xdV = -( 2^)0 5p.5(p - k) 4(5(q - k) 


S{k'^ - m^) = -— [(5(fc0 + ei^) + 5{k^ - g^)] , 

where et = v/k^ -|- m^, one can then integrate in x and y, thereafter in fcO, pO^ gO g^g obtain 


HT) = - 


(27r)3 


d^kdOpd^q 


1 




(‘5+++ -b • • • -b S -) dp^6{p - k) dq^Sici - k) 


where the S terms correspond to the 8 possible combinations of positive and negative frequency of the k, p and q 
four-momenta. Thus, we have 


5. 


+++ = + ’^T(ek)][l + nT(ep)][l +?^T(eq)] 


5_++ = (A;_,p+,g+)nT(£k)[l+ nT(ep)][l+?^T(eq)] 


jOOlOOlOO^ 

, 00 | 00 | 00 / 


-J, (fc_,p_,g_)nT(ek)?T-T(ep)»^T(eq) • 


5 ___ = Vj 

where k± = ±ek, and similarly for p and q. We can then integrate in p and q to get: 


ktiT) = - 


1 


(27r)3 


d-'k 


1 52 

8ek dp3 dqs 


1 


{S+++ H-+5-) 


(39) 


- p= —k, q=k 
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TABLE I: The integrand functions I{k,m,T) (see eq. (40)) for the correlators in (30) of a free real scalar field. 


I(k, m, T) 

smh(^){fc^e| + + 3e|(l — 4^)][cosh(^) — 1] — 

288T'^ek + 3T^[k^{l - 4^) + e|(8^ - l)][cosh(^) 

sinh“®(|i)sinh(^){fc^e| + l5T^el{l - 2f)lcosh(?A) - 1] - 5Tk 
7^ sinh“‘^(||;) sinh(^)(l — 4^) 


h 

kg 

ks 

jt 

je 

js 

A 


2Tt 


'sk sinh(^)} 


) -r <7k\^‘^ — -7^j — 1] — 2Tk £k sinti( 

i + 15 T 2 e 2 (i _ 25)[cosh(^) - 1] - sinh(^)} 


24,ek 

^ sinh-^df;) sinh(d)(8C - 1) 
sinh-^df;) sinh(d)(2C - 1) 


TABLE II: The correlators (30) calculated for a free real scalar held with vanishing chemical potential. Also shown 

the well-known expressions of p and p. 



«:(£) 

ar{x,^) 

p 

1 ^2 

30"^ 

— {rx)~^Kgirx) {rx)~^K 3 {rx) 

p 

1 ^2 

90"^ 

{rx)~^ Kgirx) 

kt 

1 .^2 1 1 t 

30^ 6 ? 

H — 1 + 24£x“^]A2(ra;) -|- 3[r(l — 8£) — 3r“^]a:“^A3(ra;)} 

kg 

90 18 ~ 3^ 

{8(1 — 5^)x~'^Kgirx) + [r(16£ — 3) — 3r“^]a:“^A3(ra;)} 

ks 

“90^ A J2 “ 2^ 

— 1 {2(1 — 2^)x~‘^Kgirx) + [r(4£ — 1) -I- A3(ra;)} 

jt 

1 (1 - 4£) 

{l-AC)X-^K2{rx) 

je 

H (8C - 1) 

|(8£ — l)a:“^A2(ra:) 

js 

(2£ -1) 

\{2i-l)x-^K2{rx) 

ly 

l35^ A Jg -I- 3^ 

1 [(12£ — 6)a:“^A2(ra;) -|- (2r -|- r“^)a:“^A3(ra:)] 


All the correlators in eq. (30) can be calculated in a similar fashion although it should be pointed out that the case 


of kt is somewhat simpler because in eq. (38) only one term in eq. (37) survived. Indeed, in general, one can have up 


to four terms associated with different sets of indices. Thus, the general correlator can be expressed as an integral of 
a sum of terms analogous to that appearing in eq. (39). 


In the massless case, T is the only energy scale and, on purely dimensional grounds, one hnds that the correlators 
are of the form For instance, integrating the eq. (39) with m = 0 one obtains: 


kt(T) = 


1 2 1 ^ 
30 6 ^ 


rj-l 4 : 


For the massive case, the integration is just a little more involved. First, the angular part of the integration in k can 
be readily carried out and one is left with expressions like: 


1 

^ y d/c I(k,m,T). 


(40) 


where the function I(k, m, T) is reported in table|T]for the various correlators. The integral in eq. (40) can be computed 
setting k = msinhy, which makes it possible to extract a factor; the integral then depends on m and T only 
through the ratio x = m/T and one is then left with an adimensional integral over y that can be turned into a series 
of type (41) involving the modified Bessel functions of the second type Kn{x): 


m 

2tt'^ 


4 oo 


'^ar{x,0, (41) 

r—1 

where, as has been mentioned, x = m/T. The final expression of functions a and k can be found in table [n| 


























TABLE III: The coefficients of the stress-energy tensor in eq. (341 calculated for a free real scalar field with 

vanishing chemical potential. 
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ar[x,G) 



^(1-60 

^ [(r^ -1- 24^a;“^)A'2(ra:) -|- 3(1 — d,G)rx~^Ki{rx)\ 

±m^T-^ + lm(l - 80 + (^ - U)T + o{T) 

f/„ 

^(1-4?) 

1(1 - A(,)x~'^K2{rx) 

(1 _20r + o(r) 

Dc 

n(6e-i) 

^ [(12 — A%G)X~'^K 2 {rx) -1- (24^ — b)rx~'^Ki{rx)\ 

m(C-|i) + (H-^)r + o(r) 



^x~‘^K2{rx) 

CT -t o{T) 

A 

1^(1-6?) 

1 [(4^ — 2 )x~^K 2 {rx) -1- (1 — AG)rx~^Ki(rx)\ 


W 


1(2^ - l)a:“^A2(ra:) 

{i-\)T + o{T) 

G 

3i(i + 6e) 

1 [(6^ — 3)x~^K2(rx) + rx~^K3(rx)] 

|m + (c-5^)r + o(r) 


With the correlators calculated, we are now in a position to write down the coefficients of eq. (32), reported in 
table III alongside with their non-relativistic limit m/T = a; ^ 1, factorized as nf{m,T) where: 


2'k'^ 


^(rx) ^K2{r 


(42) 


is the particle density at the homogeneous equilibrium. The non-relativistic limit can be extracted by simply taking 
the asymptotic expansion of the r = 1 term of each series. 


As it can be seen from the table III all the coefficients U, A, D,W,G have a finite non-relativistic limit with the 


dimension of an energy per unit volume. Consequently, as it has been mentioned, all the corrections to the stress- 
energy tensor in eq. ( |34| ) are of quantum origin as they linearly depend on h. 

The coefficient W = Aa/T^ for the massless case turns out to be in agreement with the calculation in ref. for 
^ = 0. However, unlike therein argued, we found that it has an explicit dependence on that is on the stress-energy 
tensor form. 


VI. THERMODYNAMICAL INEQUIVALENCE, FRAME DEPENDENCE AND EQUATION OF STATE 


We are now going to discuss some physical consequences of the general form of the stress-energy tensor (34) which 
we rewrite here: 


=[p + a'^Uc. + u'^u'' -[p + -h 

-l-Ad^d'^d'^ -I- + Gaui{u^j'' + 


(43) 


where the shorthands a = h\a\/cKT and uj = h\a\/KT for the adimensional scales related to acceleration and vorticity. 

The first remarkable consequence is that, as pointed out in refs. mnn], the mean stress-energy tensor in a general 
thermodynamic equilibrium depends on the fundamental stress-energy tensor operator written in terms of the quantum 
fields. This is at variance with the familiar homogeneous equilibrium, and it is made apparent by the dependence 
of the thermal functions other than p and p in table |III| on the parameter If one was able to measure one of 
the coefficients multiplying d^ or with a thermodynamics experiment, one would obtain information about the 
true, physical stress-energy tensor operator, hence on the correct gravitational theory, a conclusion already drawn in 
ref. | 21 ) . 

The second consequence is that, as it is apparent from the eq. (43), u'' = is not an eigenvector of if 7 
is non-vanishing, that is if the three vectors a,w,u (or a,oj,u) are linearly independent, as it can be seen from the 
eq. (34). This is what happens for the the rigid rotation, where a,oj and u are orthogonal to each other. In this case. 


the u vector does not coincide with the Landau definition of four-velocity, and should then be taken as defining a new 
hydrodynamical frame, dubbed the /3 frame, as it has been extensively discussed in ref. [S]. 

The third, and perhaps the most striking consequence, is that the dependence of energy density and pressure on 
the temperature and chemical potential are modified with respect to the homogeneous equilibrium case. Also, there 
are more second-order coefficients in the expansion of the stress-energy tensor than previously envisaged. Looking at 
the eq. (43) it can be realized that, with respect to the expansions presented in refs. unis], there are three new 
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coefficients, that is G, Ua, U^, and two of them imply a modification of the energy density. One could argue that they 
would disappear by going to the Landau frame. Yet, in the diagonalization of the stress-energy tensor in eq. (431, it 
can be readily shown that, retaining only quadratic terms in a and w: 


PeB = P + a^Ua + + 0{w'^) 


Peff = P + 


D, 








-I- o{w‘^), 


(44) 


where the effective pressure has been defined as the mean of the eigenvalues of the spacelike eigenvectors. Therefore, 
the energy density and the pressure coincide, in this approximation, with those in the /? frame and the coefficients 
Ua and Uw survive. One may wonder whether the modification of the energy density could be reabsorbed by a 
redefinition of the temperature other than the length of the /3 vector in the density operator in the eq. , which is 
based on the maximization of entropy with macroscopic constraints [9]. In fact, a redefinition would cure only one of 
the eigenvalues of the stress-energy tensor, unless the coefficients U, D, A, W fulfilled some preculiar relations. In all 
other cases, the relation between the eigenvalues of the stress-energy tensor, or the relation between energy density 
and pressure, in other words the equation of state Peff(Peff), is modified with respect to the homogeneous equilibrium 
case. For instance, in the non-relativistic limit of the massive case m ^ T one has, according to table |III| that the 
leading corrections are those in a^, and restoring the natural constants: 


1 mc^ _9 / 1 

Peff=sp+-—( 1 + 


24 KT' 


24 {KTf 


Peff P + 


-e-- 

3^ 8 


mc^a^n = p 


1 


2 ^ 1 \ 

3^“ Sy (KT)^ 


(45) 


where p = mn and p = nKT are the usual non-relativistic expressions for the ideal Boltzmann gas and n has the well 
known approximate expression: 


(S) 


3/2 




We note in passing that the relations (45) should hold in the case of a charged scalar field in the non-degenerate 
Boltzmann limit with a chemical potential, that is: 


fmT\ 


3/2 




and negligible anti-particle contribution. 

If it was possible to redefine T to a new T' = T + b{T)a^ such that p = mn{T') and p = T'n{T'), then the 
coefficients in the cP expansion of the functions would be the same. This can be shown by taking into account that 
dn/dT ~ {m/T^)n{T) in the non-relativistic T limit, so that 

p{T') = mn{T') mn{T) + ^(T' - T) = mn{T) (l + ^ba^) 

p(T') = T'niT') Tn{T) + n{T) (l + c. Tn{T) (l + J^ba^) 

However, it can be seen by comparing the above equation with (45) that in general this is not the case, except when 
^ = 1/4 which is neither the canonical nor the improved tensor. 

Furthermore, in general, the redefinition of a temperature would be mass dependent and it would then be trouble¬ 
some to define thermodynamic equilibrium at a common temperature of a mixture of gases. Let 

Pes{T,a,u}) = p(T'(T, a, u))), 

where p is the familiar homogeneous energy density. Expanding the new temperature in a and w the leading order 
corrections must be of the second order: 

r = T + Ta{T)a^ + TUT)oj^ + 

where Tg and Tq are proportional to the second derivatives of T'(T,a,uj) with respect to a and oj respectively. These 
unknown functions can be obtained by comparing with the equation (44): 

p + + TquP') = p + a'^Ua + 

al 
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implying 


T' = T- 




U„ 


Looking at the tables (|^ and (III), it can be realized that the coefficients of and are non-trivial functions of 
the mass and temperature. 

Going now back to the properly defined T = we observe that, in the non-relativistic limit the relation 

between the effective energy density and pressure gets modified into: 

| 2 - 


KT 


Peff - Peff- 


1 + 


3 6 


1 \ mh?\a\' 


{KTf 


Therefore, the effective equation of state depends on the acceleration besides the temperature. This could be surprising, 
but in fact in general global equilibrium all parameters, including acceleration and angular velocity play the role of 
thermodynamic variables on the same footing as temperature and chemical potential. It can be seen that in the non- 
relativistic non-degenerate limit the quantum correction to the relations ( |44[ ) and the equation of state becomes more 
important at low proper temperature, being proportional to 1/T^. Of course this applies as long as the acceleration 
is such that mh^\a\^ / {KTY ^ 1 so that the expansion method holds ^ for very low temperatures, one would have 
to take more and more terms into account and eventually the exact solution would be needed. 


VII. CONCLUSIONS 

In conclusion, we have demonstrated that the relativistic stress-energy tensor in general states of global thermo¬ 
dynamic equilibrium features quantum corrections with respect to its ideal form ([^ depending on the local values of 
acceleration and vorticity, besides proper temperature and chemical potential. We have calculated the coefficients of 
the additional terms of the stress-energy tensor in the appropriate quantum statistical framework at the second order 
of an expansion in the parameters Ha/cKT and Tloj/KT for the simplest case of a real scalar field. We have found 
that more terms exist with respect to previous assessments; our calculated coefficient W for the real scalar field agrees 
with previous calculations [6]. 

We have emphasized three major physical consequences of this finding: 

1. The coefficients explicitely depend on the form of the quantum stress-energy tensor operator, what was already 
argued in refs. [Hill]. 

2. The effective energy density - defined as the eigenvalue of the stress-energy tensor - is also modified by terms 
involving acceleration and vorticity which cannot be reabsorbed by means of a redefinition of the temperature. 

3. The equation of state and the relation between effective pressure and effective energy density are also modified 
by the presence of vorticity and acceleration. 

In principle, these findings could be extended to matter in local thermodynamic equilibrium in flat spacetime, as 
well as to matter in global/local equilibrium in a curved spacetime. In this case, it is well known that ^ in eq. (§ 
must be a Killing vector which can have a non-vanishing exterior derivative 9^/3,^ — and, consequently, additional 
terms of the stress-energy tensor with respect to its ideal form ([^. This might be of phenomenological relevance for 
the study of the equilibrium of self-gravitating objects. 
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Appendix A: TZ expansion 


To derive the espressions of , we can disregard, for the sake of simplicity, the conserved charge in (221 for it 
commutes with both P and J operators. Defining 

A = B = , 

and applying the known Poincare algebra relations, we find: 


(Al) 


[[A,B],A]=0, 

[[A,B],B]=l3^^vj^,u;p,g'^pp‘^, 

[[[[A, B],B ],...], B] = -(i)”/3^(ro -zu- ...■ zu)^^P''. 

Now, using the Baker-Campbell-Hausdorff formula to expand exp[A + B] and retaining only the non-vanishing terms; 
taking into account that any commutator involving A, B or commutators thereof, will in turn commute with A, being 
proportional to four-momentum operators, we obtain: 


ni3,w) = e^+^ cs e^ e^e^I^’^1 ... 




and its expansion up to second order in B (which is tantamount to a second order in zu) reads: 

P(/3, zu)^e^+{B+ ^[A, B])e^ + ^ Jb[A, B] + ^[A, B]B + ^[A, B]^) e 


6 ' 


8 ' 


(A2) 


(A3) 


where advantage has been taken of the fact that exp[A] commutes with both the commutators [A, B] and [[A, B], B]. 
Now, by using the relation: 


e-^Be^ = B - [A,B] 


which is a known corollary of the Baker-Campbell-Hausdorff formula for our case, the eq. (A31 can be rewritten as: 

A 


1 


1 


1 , 


1 , 


P(/3, tn) ~ e^ + e^(B - -[A, B]) + e^ ( -B^ - -B[A, B] - -[A, B]B + -[A, B] 


(A4) 
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We can now take the half-sum of (A31 and (A41 to obtain: 




(A5) 


putting the expressions of A and B in eq. (Al) in the eq. ( |A5[ ) one can read off the operators in eq. (23), which are 
quoted in eq. (1^. 


Appendix B: Calculation of angular momentum-stress energy tensor correlators 

The density operator ([^, which is used to calculate the mean values denoted as {)p can be written as 
where A is the Lorentz transformation turning /3o = (1/T, 0 ) into /3. Hence, the mean value of a general tensor can 
be expanded as: 

= (bi) 

where ( )t, as has been mentioned in the text, stands for the mean value with the density operator ^ exp[—/3o • P] = 
^ exp[—H/T], Note that: 










(B2) 


Since exp[—H/T] is invariant by rotation, only scalars under spatial rotation, either components or contractions of 
the tensor may have a non vanishing value. Furthermore, we assume that the hamiltonian operator is 

symmetric under parity and time reversal transformations, so that also pseudoscalars and scalars which are odd under 
time reversal will have vanishing mean value. 

For instance, for a symmetric tensor operator one can write: 


= Si;S^A + A^^'B 


(B3) 


where is the transverse projector in the rest frame, i.e. A^‘' = g^" — SqSq. 
well known: 


Of course, the eq. (B3) becomes the 




by usi ng th e (Bl) and ( |B2[ ). The coefficients A and B can be calculated from the mean values selecting the components 
in eq. (B3| which make all terms vanishing except one. Thereby: 


A = 


B = 


This general procedure can be applied to the calculations of tensor s of any rank. Indeed, in view of eq. (B2|, 
anytime a time component or 0 index is selected in in eq. (Bl| a Sq will appear eventually turning into 

a u after boosting, while for a space contraction of indices a A projector will, eventually turning into a A like in 
eq. (IbI. 


We can first appl y th e above argument to the calculation of 0 . By using the decomposition (27) and 

taking into account (Bl), the only possible non-vanishing contributions read: 




Yet, they both vanish because they are odd under time reversal and parity respectively. No scalar can be formed with 
{ and so the mean value of the angular momentum vanishes too. 

Let us now move to the more complicated case of correlators involving two angular momentum operators, starting 
from: 
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In the rest frame, K'^} is a symmetric tensor with vanishing time components, so it has one spin-0 component 
obtained with the contraction of the indices p and a and one spin-2 component under rotation which can be obtained 
by applying the projector: 


^ ^ 

to the tensor itself. In order to construct a rotation singlet, we need to combine the above components with the 
corresponding components of . The spin 0 components can only contract with and its spatial trace, so one 
obtains two contributions: 


whereas the contraction of the spin 2 component of with the one of T gives rise to: 




“^Red.^^, k'^y, fP'')T = -AP'^S^S^kt + AP'^AP'^ke + ( Ap^A'^'' + AP''A'^p - ‘^A^'^AP'^ ) (B4) 


Altogether 


which in the observer frame reads: 

^Red^'’, k'^}; fP‘')p = -AP’^uPu''kt + AP’^A^^'ko + (^Ap^A'^'' + AP''A'^p - ? ks 

To find a compact expression of the coefficients kt,ke, kg one can select the indices making all terms on the right hand 
side of (B4) vanishing except the one of interest. One can check that all indices in the definitions (30) are properly 
chosen (notice how in kg we avoided the symmetrization in /i -f-)- associated with the anti-commutator since we know 
that the anti-symmetric part will not contribute). For the jt,je,js the procedure is precisely the same outlined above 
with the replacement K ^ J. 

In fact, the correlator ({AT^, k}', TP'')t is a somewhat special case because {KP, J'^} is odd under parity and time 
reversal. Therefore, the only non-vanishing contraction is between the two spin-I components of the tensors {AT^, 
and TP'' respectively. The spin-1 components can be obtained by means of the projectors 

1 


pP'T — —ZeOpCTT 

(1/3 “ 2^ S0a/3T 


and 


Pkp = i^5PA''Jl + 5fAPJl 

respectively. Hence: 

RediFd -h d = -b SIc’^p’^pS^) 


which, once boosted, reads: 


RediFd ry TP'')p = {u^e^P^'^uP + 






